###########################CARREGAR PACOTES##############################
#########################################################################
require(RcmdrMisc) #pacote para carregar arquivos do Excel
require(agricolae)
require(dae)
require(ExpDes)


#Mudar diretório
setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
source("funcoes.txt")



#########################################################################
############################Análise do DQL###############################
#########################################################################

#########################################################################
############################Análise do DQL###############################
#########################################################################

#Mudar diretório
setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")

source("funcoes.txt")


##Importar dados
dados=readXL("anestesico.xls",sheet = 1,stringsAsFactors=TRUE)
attach(dados)
str(dados)

attach(dados)
str(dados)
linha=as.factor(linha)
tapply(trat,list(linha,coluna),sort)

######Análise de Variância#####################################3
model=aov(vr~coluna+linha+trat)
anova(model)
cv.model(model)

#####PRESUPOSIÇÕES DA ANOVA
#Gerar os erros padronizados do modelo
e=rstandard(model)



#Teste de Shapiro-Wilk - NORMALIDADE
shapiro.test(e) 

#Teste de Oneill-Mathews - Homogeneidade de Variância
VR=vr
Linha=linha
Coluna=coluna
Trat=trat
oneillmathewsDQL(VR,Trat,Linha,Coluna)

##Teste de Durbin-Watson - Independencia
set.seed(12345)
durbinWatsonTest(model)


###Teste de aditividade dos efeitos de tratamento, linha e coluna
model1=aov(VR~Coluna+Linha+Trat + Error(Linha*Coluna))
tukey.1df(model1, dados, error.term= "Linha:Coluna")



####Anova com Teste de Tukey####################
latsd(Trat, Linha,Coluna,VR, quali = TRUE, mcomp = "tukey", sigT = 0.05, sigF = 0.05)

####Anova com Teste de Scott-Knott####################
latsd(Trat, Linha,Coluna,VR, quali = TRUE, mcomp = "sk", sigT = 0.05, sigF = 0.05)